Two-Point Green's Function in VT- Symmetric Theories 
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Abstract 

The Hamiltonian H = ^p 2 + ^m 2 x 2 +gx 2 (ix) s with 5, g > is non-Hermitian, 
but the energy levels are real and positive as a consequence of VT symme- 



^ ■ try. The quantum mechanical theory described by H is treated as a one- 

dimensional Euclidean quantum field theory. The two-point Green's function 

00 

for this theory is investigated using perturbative and numerical techniques. 
The Kallen-Lehmann representation for the Green's function is constructed, 
and it is shown that by virtue of VT symmetry the Green's function is entirely 
D . real. While the wave-function renormalization constant Z cannot be inter- 



preted as a conventional probability, it still obeys a normalization determined 
by the commutation relations of the field. This provides strong evidence that 
the eigenfunctions of the Hamiltonian are complete. 

In this Letter we study the quantum mechanical theory described by the Hamiltonian 

H = ip 2 + imV + gx 2 {ix) s {5, g>Q). (1) 

We treat this theory as a one-dimensional self-interacting scalar Euclidean quantum field 
theory whose Hamiltonian is 

H = i(V0) 2 + im 2 2 + g^ 2 ^ (5, g > 0) (2) 

and study the two-point Green's function for this theory. 



The quantum mechanical theory in (|l|) has been studied in great detail in the past few 
years because, while it is not Hermitian, it does possess VT symmetry. (Under parity V, 
x — > —x and p — ► — p; under time reversal T, x — > x, p —>■ —p, and i — > —i. Thus, the 
operator VT commutes with the Hamiltonian H.) It appears that, as a consequence, all of 
the energy levels are real and positive PHI^C- The first rigorous proof of the reality of the 



spectrum for the massless case m = is given in Ref. [|13|]. Direct numerical evidence for the 
reality and positivity of the spectrum can be found by performing a Runge-Kutta integration 
of the associated Schrodinger equation [|lj . Alternatively, the large-energy eigenvalues of the 
spectrum can be calculated with great accuracy by using WKB techniques [p],^|,|15|. One 
can also demonstrate the reality and positivity of the spectrum by calculating exactly the 
spectral zeta function (the sum of the reciprocals of the eigenvalues). This was done by 
Mezincescu JFj and Bender and Wang 0. 

Non-Hermitian quantum mechanical PT-symmetric theories have been studied exten- 
sively, but not much work has been done on the corresponding scalar quantum field theories. 
Three major studies have been published. The first [[ll]] shows how to solve the Schwinger- 
Dyson equations for such theories. The interesting result is that the beta functions for the 
ig<f) 3 and the —g4> 4 theories have the opposite signs from those of the conventional g<p 3 and 
g<j) 4 theories. Thus, in four dimensional space these non-Hermitian theories are asymptot- 



ically free and are nontrivial. The second paper |L7| is a nonperturbative calculation in 
D- dimensional space of the one-point Green's function G\ in such theories. It is quite re- 
markable that G\ = (O|0(O)|O) is nonvanishing for all 5 > 0. This is true even for the case 
6 = 2, which corresponds to a — g4> 4 field theory. Hence, this model might possibly be useful 
to describe the physics of the Higgs sector. The third paper |18[ concerns bound states 



in such theories. It is shown that for sufficiently small coupling constant, a —g4> 4 theory 
possesses n-particle bound states for any n. In fact, there is a range of g for which there 
are just two-particle and three-particle bound states but not bound states of four or more 
particles, which is somewhat reminiscent of the way in which quarks form bound states. 
In this Letter we consider the field theoretic connected two-point Green's function, which 



in Minkowski space is defined in general by 

G Minkowski (;E _ y) = ( O |T[0(a:)0( 2 /)]|O) connected , (3) 

where T represents time ordering. To discuss this equation in the context of quantum 
mechanics, which is equivalent to quantum field theory in one time and no space dimensions, 
we incorporate the time ordering by using the Heaviside step function 9{x): 

(O|T[0(x)0(y)]|O> = (0|<Kx)<K2/)|0)#(x -y) + (O\<l>{y)<Kx)\O)0(y - x). (4) 

We evaluate (0\<p(x)<p(y)\0) by inserting a complete set of states: 

oo 

(0\<j>{x)<f>(y)\0) = £<O|0(a:)|7i)H^/)|O) 

n=0 

oo 

= E (0\e~ iHx (p(0)e iHx \n) (n\ e - iHy (f)(0) e iHy \0) 

n=0 

oo 

= £[(O|0(O)|n)]V (2 ^-- E * ))(a, - 1 ' ) . (5) 

n=0 

Thus, if we include the effect of time ordering, we obtain 

oo 

(0|Tfe(aMl/)]|0) = Y,mm\n)} 2 e^- E ^-y\. (6) 

n=0 

In this sum the n = term [(O|0(O)|O)] 2 , which is the square of the one-point Green's 
function, represents the disconnected part of the two-point Green's function. Thus, the 
connected two-point Minkowski space Green's function is given by 

oo 

G Mi„kowski( x _ y j = ^[( | ( /,(0)|n)] 2 e i(B "- So)|1 ^ 1 . (7) 

n=l 

In Euclidean space we replace the i in the exponent of (J7h by —1 and get 

oo 

G 2 (x-y) = E[(O|0(O)|n)] 2 e-^-^)l^l. (8) 



n=l 



We perform a Fourier transform to find the Green's function in momentum space: 



where we use the identity g2 2 ^ 2 = /f^ dxe~\ x ^e %Ex . 

The difference of the first-excited-state energy E\ and the ground-state energy E is 
defined as the renormalized mass M. It is traditional to define the wave-function renormal- 
ization constant Z as the coefficient of the one-particle pole E 2 +M ? ■ In quantum mechanics 
we have an infinite sum of poles instead of a cut. Thus, we define the nth wave- function 
renormalization constant Z n by 

Z n = 2(E n -E )[(0\<j ) (0)\n}} 2 . (10) 

Finally, we express the matrix elements above in terms of the usual coordinate-space energy 
eigenfunctions ip n (x), which satisfy the time- independent Schrodinger equation 

FVnO) = E n tfj n {x). (11) 

We thus obtain a representation for Z n as a ratio of integrals: 

? = 2(£ n - E Q ) [J dx x^ (x)^ n (x)] 2 
n [Jdx^(x)}[Jdx^ n (x)} ■ [ ' 

We have not seen this general result before in the literature. 

We now make a crucial observation. For ordinary real Hamiltonians the solution to the 

Schrodinger equation ([0]) is real, so it is not necessary in (0) to use complex conjugates 

of wave functions. That is, we have written ipn( x ) an d not \i/j n (x)\ 2 . In our upcoming 

discussion of the complex Hamiltonian ([I]) we will continue to use the formula (12). It 



would be wrong to use complex conjugates of wave functions because for non-Hermitian 
'PT-symmetric Hamiltonians the statement of orthogonality for the wave functions is[| 

dxi[) m (x)i/) n (x) = (-l) n S mn . (13) 



It is important that the integrand in (13) not contain the complex conjugate of ip m {x) because 



it is not an analytic function of x. The integrand of fll3| ) must be analytic because we will need to 
have the freedom to move the contour of integration in the complex-x plane as the parameter 5 in 
the Hamiltonian (|l|) varies. Thus, we do not distinguish between the states \n) and (n\. 



In terms of Z n the expression for the momentum-space two-point Green's function is 

oo y 

where M n = E n —Eo. This formula is the quantum mechanical analog of the Kallen-Lehmann 



representation [1^] for the two-point Green's function. 



The wave-function renormalization constants Z n must satisfy the general constraint 

oo 

5> B = 1. (15) 

71=1 

To derive this result we evaluate the sum by substituting fllCf): 

oo oo 

E^ = E 2 (K-^o)[(O|0(O)|n)] 2 
E2(^- J E o )[(O|0(O)|n)] 2 



n=l n=l 

oo 



n=0 
oo 



E [<O|0(O)tf |n)H0(O)|O> - (O|#0(O)|n)<r#(O) |0) 



n=0 



+ (0 |0(O)|n)(n|#0(O)|O> - (O|0(O) |n) <n|^(0)iJ |0)] 
= (0| [0(0)^0(0) - ^(0)0(0) + 0(O)#0(O) - 0(O)0(O)tf] |0), (16) 

where we have assumed the completeness property 1 = J2^=o \ n ){ n \ f° r the states |n). We 



then reduce (|T6|) using the equal-time commutation relation of the field 0: 



E^n = (O|[[0(O),tf],0(O)]|O> = (O|#(O),0(O)]|O> = (0|z(-z)|0) = (0|0) = 1. (17) 

n=l 

We can calculate the two-point Green's function using several different kinds of per- 
turbative methods. For example, we can sum the one-particle irreducible Feynman graphs 
having two external legs. Alternatively, we can expand the wave functions and energies in 
Rayleigh-Schrodinger perturbation series and then use (|12|) and fll4|) . For a one-dimensional 
conventional g<j) 4 theory whose Hamiltonian is 

H=^p 2 + i m V + gx 4 {g > 0), (18) 

the Euclidean-space Feynman rules are —24g for a vertex and E i+ m $ for a line. The results 
for the first few wave function renormalization constants are 



^ 1 9 2 297 3 138405 4 1015137 5 991058565 6 nnAnnroAcn nrr 7 

Z 1 = 1 e H e d e 4 H e 5 — e 6 + 79420584690 + ver256e 7 + 

8 8 128 32 1024 

9 , 297 o 138105 , 2015199 = 974210895 fi 30738579111 7 

Z 3 = -e 2 e 3 H e 4 e° H e 6 e 7 H , 

d 8 8 128 64 1024 1024 



Z5 "32 e ~~6^ e + ^12 _ C 1024^ C + '"' 

2205 fi 222705 7 

Z 7 = e b e 7 H , 

7 256 128 

where e = g/m 3 . For this theory the ij2 n vanish by parity. Evidently, the constraint ([IB]) is 
satisfied perturbatively: 

Z 1 + Z 3 + Z 5 + Z 7 = l + 0(e 8 ). (20) 

What happens for the non-Hermitian — g^ 4 theory that is obtained when we set 5 = 2 
in ([!])? A perturbative calculation in powers of g gives exactly the same results as above 
except that e must be replaced by — e in (|i~9|) . Thus, the perturbation expansions now do 
not alternate in sign. As a result, one might think that each of the Z 2n +i would have an 
imaginary part because the series in ( |i~9D are divergent. (That is, one might expect there to 
be a cut on the positive axis in the complex-e plane.) However, this is not the case. Indeed, 
it follows easily from ([12]) that Z n is real for all n. The key point is that the wave functions 
are all eigenstates of VT ||, and thus exhibit the symmetry [^(x)]* = ip(—x). 

Operationally, the reason that Z n is real is that there is an exponentially small nonpertur- 
bative contribution that exactly cancels the discontinuity across the cut. This contribution 
occurs in all the Z n , so it is no longer true that Z 2n = 0.0 

Using Runge-Kutta we have found the numerical solution to the Schrodinger equation 
(|TT|) and computed the first eight Z n for a — g<p A quantum mechanical theory using (|l~2|). We 



2 One must be careful to avoid drawing wrong conclusions on the basis of perturbative calculations. 
For example, a perturbative calculation of the one-point Green's function G\ gives to all orders 
in powers of the coupling constant g. However, the correct result is nonzero, pure imaginary, and 
exponentially small |17|]: G\ ~ — iKm 3 / 2 g~ 2 ^ 3 exp[— m 3 /(3<7)] as g — > + , where K is a positive 
constant. This nonperturbative result is due to a soliton solution of the classical field equations. 

6 



find that the Z n alternate in sign with Z\ being positive. When e = g/m 3 = 0.01, we can see 
the nonperturbative behavior of Z 2n ; these numbers are nonzero and exponentially small: 



Zi 


0.9998757421 




z 2 


-0.2341326909 x 10" 


-17 


z 3 


0.1649676955 x 10" 


-3 


z± 


-0.6881796661 x 10" 


-15 


z 5 


0.8686216150 x 10" 


-7 


z 6 


-0.8216102337 x 10" 


-12 


z 7 


0.2227856879 x 10" 


-9 


z» 


-0.1027951293 x 10" 


-10 



(21) 



Apparently, the Z 2n increase until they reach a maximum and then decay to zero. The sum 
of the renormalization constants is unity in accordance with fllS]) . 

When e = 0.1, the partial sums of the Z n approach 1 in an oscillatory fashion: 



n 


Z n 


ELi z k 


1 


1.1793780652 


1.1793780652 


2 


-0.2117900432 


0.9675880220 


3 


0.0362453240 


1.0038333461 


4 


-0.0042084418 


0.9996249043 


5 


0.0004074167 


1.0000323210 


6 


-0.0000348628 


0.9999974582 


7 


0.0000027315 


1.0000001897 


8 


-0.0000001996 


0.9999999901 



(22) 



It is difficult to prove completeness when the Hamiltonian is not Hermitian, but the fact 
that the numerical sum of the renormalization constants is unity is powerful evidence that 
the states are complete. This is the first time that direct evidence for the completeness of 
the states of H in ([!]) has been presented.^] 

3 Indirect evidence for completeness in the massless (m = 0) case was found in Ref. || . 



7 



In conclusion we point out that many of the results we have found here for one- 
dimensional Euclidean space immediately extend to the theory characterized by the Hamil- 
tonian density (|2|) in D-dimensional Euclidean space. In particular, we can show that the 
2n-point Green's functions are all real and that the (2n + l)-point Green's functions are 
all pure imaginary. This result follows from the VT symmetry of the path-integral repre- 
sentation of the Green's functions. A detailed study of the wave function renormalization 
constant for D > 1 will be discussed in a future paper. 
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